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Potential-density pairs 
and vertical tilt of the stellar velocity ellipsoid 

O. Bienayme 



0^ 
O 
O 
(N 

< 



< 

o 

Oh 
I 

o 



> 

O 



o 

On 
O 



X 



Universite de Strasbourg, CNRS, Observatoire Astronomique, France 
09/03/2009 

ABSTRACT 

We define new potential-density pairs and examine the impact of the potential flattening on the vertical velocity ellipsoid tilt, 6. By 
means of numerical integrations and analytical calculations, we estimate 5 in a variety of galactic potentials. We show that at 1 kpc 
above the Galactic plane at the solar radius, d can differ by 5 degrees, depending on whether the dark matter halo is flat or spherical. 
This result excludes the possibility of an extremely flattened Galactic dark halo. 
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1. Introduction 

The velocity ellipsoid shape can be determined from the stellar 
kinematics in our Galaxy, see for instance Chereul et al. d 19981 1 
from Hipparcos data, Soubiran et al. (2008 ) from distant red gi- 
ants or Siebert et al. (2008 ) from Rave data. It is known that 
the vertical tilt of the velocity ellipsoid is related to the galac- 
tic gravitational potential and that, in the peculiar case of the 
Stackel potentials, it only depends on the gravitational poten- 
tial. Probably because of the lack of data, very few works have 
examined this question since the pioneering work of Ollongren 
(11962b . We can note the work by Kuijken and Gilmore (11989b 
that proposes an approximate determination of the inclination at 
the solar position above the galactic plane and that points out the 
necessity to consider this inclination to determine the K, force 
out of the galactic plane. There are also the analytical studies 
by Cuddeford and Amendt (199n of the Jeans equations to the 
fourth order and their analytical approximation of the vertical 
tilt. Many other studies concern Stackel potentials for which the 
ellipsoid orientation is easily obtained. 

Thus for realistic galactic potentials and from orbit integra- 
tion, both Kuijken & Gilmore (1989) and Binney & Spergel 
( 11983b show that at 1 kpc out of the galactic plane above the 
Sun, the velocity ellipsoid points towards a point located at 5 to 
10 kpc behind the galactic center (see also Kent & de Zeeuw, 
[1991 1 and Shapiro et al.. [2003T l. 

In Section 2, we define a new potential-density pair to mea- 
sure (Section 3) the vertical tilt of the velocity ellipsoid numer- 
ically within a spheroidal potential by varying its flattening. In 
Section 4, we generalize the potential-density pair to a three- 
parameter family with properties similar to the Miyamoto and 
Nagai potentials, but with flat circular velocity curves at large 
radii. In Section 5, we define a simple expression to estimate the 
tilt and examine its accuracy. 



2. A new potential-density pair 

We make the definition of the Mestel disk more general by uti- 
lizing the scale-free axisymmetric potential 



(1) 



It generates, within the plane of symmetry z=0, a constant cir- 
cular velocity curve of amplitude v^. The resulting density is al- 
ways positive: 



p(R,z) 



47tG (qlRl + ^2) 



1/2 



2^3/2 ' 



(2) 



Figures T]|2 show the isodensity and isopotential contours for 
some values of q. The axis ratio of isodensity contours is q^^^ 
and for the isopotentials it is (1 + q)/2, never less than 1/2. It 
can be shown (Toomre, 1982 ) that any scale-free potential at the 
spherical limit (obtained here with q=l) is the potential of the 
singular isothermal sphere, while the flat disk limit (^=0) is the 
Mestel flat disk. With q>l, the Eq. [T] potential is prolate and 
at the limit q=oo the potential and density distributions do not 
depend on the z coordinate (cylindrical distribution). 
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Fig. 1. Contours of isopotentials with ^=0.8 and 0.2. 
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Fig. 2. Contours of isodensities with ^=0.8 and 0.2. 



3. Velocity ellipsoid tilt 

Since the Eq.[T] potential is scale-free we simply need to deter- 
mine the velocity ellipsoid orientation versus z at some given 
radius R. We build the velocity ellipsoid numerically from orbit 
integrations (2. 10^ orbits with ct^Io-r ~1/2 at z=0). The veloc- 
ity ellipsoid orientation for these distribution functions is shown 
(Fig. [3]l, for R-l, versus z within the meridional plane for po- 
tentials with flattenings q=Q.3, 0.7. Other approaches could be 
tried: see for instance, the analysis by Evans et al. (il997 ) from 
the Jeans equations for scale-free potentials with a flat rotation 
curve by applying an approximate third integral (see also for 
scale-free potentials: Qian et al. 1995, Evans & de Zeeuw [r992l 
Qian|T992l Hunter et al. [19841 de Zeeuw et al. [T996l Evans et al. 
fT997T l. 

Figure |4] shows the variation in the tilt angle 6 at R-E and 
z= 1 above the mid-plane versus the coeflicient q. 6 varies by 5 
degrees from 2.2 to 7.1 degrees depending on whether the den- 
sity distribution is flat or spherical. This may be compared with 
the recent determination by Siebert et al (2()08) of the inclina- 
tion of the velocity ellipsoid at 1 kpc below the Galactic plane. 
From RAVE data (Zwitter et al. 2008 ), they find an inclination 
6 = 7.3 ±1.8 degrees towards the Galactic plane. They show it 
implies a moderate flattening of the dark halo by exploring the 
velocity ellipsoid tilt within various galactic potentials (Dehnen 
& Binney |1998! l. Within the potential of the Besangon galac- 
tic model (Robin et al 120031 (probably the best current Galactic 
baryonic mass model) that includes a spherical dark halo, we 
also determine the velocity ellispoid tilt from numerical orbit in- 
tegration and obtain 6 ~ 6 degrees (Bienayme et al. 12009 1) . To 
summarize, within a spherical Galactic potential, the tilt at 1 kpc 
above the plane at the solar radius is about 7 degrees, and an ex- 
treme flattening of the dark matter component is excluded, since 
the tilt would be close to 2 or 3 degrees (Fig. [4]). 




Fig. 3. Vertical tilt of the velocity ellipsoid 6 in degrees versus 
z at R - 1 for models with q = 0.3 (left) and q - 0.7 (right). 
Black line: analytical estimate (see Section 5), red line: result 
of numerical orbit integration, dotted line: tilt within a spherical 
potential. 



Fig. 4. Velocity ellipsoid inclination in degrees versus potential 
flattening q at /?=8 and z- 1 . Dotted line: the expected inclination 
from Section 5, crosses: numerical determinations from orbit in- 
tegration. 



4. Other potential-density pairs 

The Mestel disk potential can be generalized by introducing a 
core radius a. It gives (Evans & Collett [1993] and quoted refer- 
ences Rybicki, 1974, and Zang |1976] l 



(piR, z) = In ( y/R^ + (a + \z\)^ -h | a -h |z| |) , 
with the resulting surface density 



and the rotation curve v^., flat at large R, is given by 

fi>2 

VciRf 



We generalize these Mestel disks with a core by applying the 
Miyamoto and Nagai ( 1975) method to generalize the Plummer 
spheroid and Kuzmin disk (see also Baes .2009j . Thus we define 
the potential 



<f>(R,z)^\n 



R^ + + y/b^ + z^j + j a h- yjb^ + ^2 



(3) 



It can be verified that this potential generates a positive 
spheroidal density distribution. With b-Q it may generate both 
a flat disk and a spheroid. The sum (a + b) constrains the size of 
the core, and b the thickness of the disky part of the spheroid. 

Potentials given by Eqs.[T]and[3]are peculiar cases (^ = or 
a - b - Q)of the more general potential 



0(^,z) = 
In 



^R^ + [a+ yjb^ + z^f + ^q^R^ + {a+ ^b^ + z^f^ . 



(4) 



and the corresponding density is (if b ^ 0) 

q^m + b^n 



{R^ + a^f^ {q^R^ + 0-2) 



2^3/2 



(5) 
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Fig. 5. Isolevels for the density distribution resulting from the potential given by Eq.|4]with (a,b,q)-{5, 2, 1.3), (5, 2, 0.7), (0, 5, 
0.5) from left to right. 



with 



n = {r^ +13^ +I3a)a^ 



a — a + 



VP 



(6) 



If a = 0, it may be shortly written as 



I q^(R^ +b^+ z^f + [{q^ + l)R^ + 2b^ + 2z^) b^ 



(R^ +b^+ z2)3/2 (^2/^2 + 1,2 + j;2) 



2^3/2 



.(7) 



When ^7 = and ^ 9^ 0, it is the combination of a spheroid and a 
flat disk 



1/2 



1 q^{R^ + (a + \z\f) 



1 



IttG (R^ + a^'Hq^R^ + a^) 



1/2- 



(8) 



(9) 



With q = 1, this is the potential-density pair defined by Brada & 
Milgrom ( fT995t . 

More generally, we obtain a spheroid potential producing a 
rising velocity curve, flat at large radii, the core radius depending 
only on q and on a + b. The parameter q modifies the flattening 
at large R and b the thickness of the disky part of the density 
distribution. 

Thus, this generalization introduces a new family of 
potential-density pairs with general properties similar to the 
Miyamoto and Nagai potentials, but with flat circular veloc- 
ity curves at large radii. Figures |5 illustrate isodensity contours 
for three such potentials. To conclude, we recall another three- 
parameter family of potential-density pairs obtained by Zhao 
(fT996b . 

5. Anaytical estimate of the vertical velocity 
ellipsoid tilt 

We examine the reliability of the Cuddeford & Amendt formulae 
(Eqs. 90-91 and ElO, 19^^ and Amendt & Cuddeford (Eq. 104, 
[1991) that predict the vertical tilt of the velocity ellipsoid within 
an axisymmetric potential, and we propose another formula that 
is more accurate at greater distances from the plane of symmetry 



z = 0. For various potentials, we compare the expected tilt from 
these two expressions with the tilt obtained from numerical orbit 
integrations. 

Cuddeford & Amendt ( |1991| l analyze the consecutive mo- 
ments of the Boltzmann equation up to the 4th order by expand- 
ing these moments in Taylor series assuming that cr,,/Vcircuiar is 
small. Combining these moment equations, they derive expres- 
sions for the velocity moments and obtain an approximate ex- 
pression for the tilt angle, tan 25 - , which only depends 

on the gravitational potential. They show that their expression 
is exact in the case of axisymmetric Stackel potentials, peculiar 
potentials for which the tilt only depends on the potentials but 
not on the exact distribution function. According to their initial 
hypothesis, the formula must be generally valid when the con- 
sidered distribution of stellar orbits covers a sufficiently small 
domain (expecting that the covered domain can be approximated 
with a Stackel potential). This should be the case of nearly cir- 
cular orbits with small vertical or radial oscillations. 

The simplicity of the result obtained by Cuddeford & 
Amendt (119911) makes it extremely attractive. For instance 
Vallenari et al. (2006) use it to improve their recent galactic 
model of star counts and kinematics, which is dynamically self- 
consistent locally. 

Inspired by their work, we propose a more direct expression 
derivated from the generic equation (Ollongren 119621 Eq. 7.2) 
that defines axisymmetric Stackel potentials: 



2 2\ , {<PR,R - 4>z,z) + 3 (Z0« - R<P:) 

z]+ , (10) 



with a positive sign for prolate spheroidal coordinates, and a neg- 
ative one for oblate ones (see de Zeeuw. [T9851 for a detailed de- 
scription). 

For a given potential 0, we substitute zo from Eq. [TOjin the 
expression for the velocity ellipsoid tilt 6 given by Hori & Liu 
(119631 1. which is exact in the case of Stackel potentials: 



tan 26 ■ 



2Rz 



R^- 



(11) 



Combining Eqs. [T0| and [TT| gives an estimate of the velocity 
ellipsoid inclination within axisymmetric potentials. We note 
that the Amendt & Cuddeford (1991) formula could have been 
obtained in the same way, just by differentiating the Ollongren 
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equation with respect to z and replacing it within Eq.[TT] 

A noticeable difference between their expression and ours is 
that their formula is defined at z = and is exact for the Stackel 
potential only at z = 0, while our expression remains exact (for 
Stackel potentials) at any z- This may be why our formula is in 
better agreement with the results of the numerical explorations 
described below. 

Now, to verify the reliability of the estimated tilt angle 5 
given by Eqs. [TOpT and the one given by Amendt & Cuddeford 
(11991b formulae, we measure the tilt numerically within a series 
of potentials. For this purpose, we build stationary distribution 
functions with a library of 2 10^ orbits for each model using a 
7th-order Runge-Kutta (Fehlberg |1968l ) with a relative accuracy 
e - 10"'^. Each orbit is integrated for 80 rotations, the initial 
conditions being drawn from a Shu distribution function (we also 
tried a Dehnen disk distribution function) with erg / Vc,>. =0.4 and 
o'zl^cir. - 0.2 (we also consider the 0.2-0.1 and 0.6-0.3 pairs). 
One point is selected per orbit in the last 40 rotations. 

Figures |6] show the inclination 6 versus the distance z from 
the plane of symmetry in the case of the Brada and Milgrom 
( 11995b disk-i-halo potential with a core radius a = 3. The tilt S 
is plotted versus z at four different galactic radii R-\.5, 4.5, 7.5, 
and 10.5. The dark line is the analytical estimate from Eqs 
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1 1 while the red line is the result of numerical orbit integrations. 
The agreement between numerical measures and the analytical 
estimate is satisfying in the range of considered positions R and 
z. The estimated vertical tilt is also shown in Figure|4] at position 
R - % and z - 1, <5 is plotted versus the flattening q of the Eq. [T] 
potential. The agreement is satisfying for q between and 1 . 




1 




Fig. 6. Velocity ellipsoid tilt angles in degrees versus z within the 
Brada & Milgrom (11995b Disk-i-Halo potential at four galactic 
radii R-1.5, 4.5,7.5 and 10.5 (from top left to bottom right). Red 
irregular line: numerical determination. Dark continuous line: 
analytical estimate. Dotted line: tilt within a spherical potential. 



We explored other potentials: a Kuzmin disk (a Stackel po- 
tential to test our numerical process), the Mestel disk, the Brada 
and Milgrom potential, a logarithmic flattened potential, com- 
bined Kuzmin-Hlogarithmic halo potentials. For these potentials 
our analytical estimate give excellent results in agreement with 
numerical determinations at a precision better than 1 degree at 
zIR = 0.15. 

With the slightly flattened spherical logarithmic potentials, 
the Amendt & Cuddeford (1991) formula is also an excellent 



estimate for the tilt (1 percent difi'erence at zIR - 0.15 with a 
potential axis ratio ^=0.9). However, with the Mestel potential, 
it fails at any z (a null tilt is predicted), while it works correctly 
with the Brada and Milgrom potential, but only at radius close 
to or smaller than the core radius. 

In conclusion, our analytical estimate of the vertical tilt of 
the stellar velocity ellipsoid is accurate for various potentials 
and vertical distances less than 0.15/?. It also shows that at 1 
kpc above the Galactic plane at the solar radius, 5 varies by 5 
degrees depending on whether the potential is spherical or the 
dark matter component is flat. 

However, for an exponential disk, this estimate fails in the 
range < 2 or > 4, in units of the scale length. The situation is 
similar for the potential given by Eq.[T]for "flattening" of ^ = 2 
and neighboring values. An explanation for this deficiency is that 
the denominator of the r.h.s. of Eq. [TTjis zero. In that case, the 
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resulting value for zo varies strongly with R or z- Since Eq 
not a Stackel fit, its domain of validity is determined by the con 
dition that the estimate for zo from Eq. 1 1 does not vary strongly 
over the considered (R, z)-domain. For deviations from this limit, 
more sophisticated methods of fitting potentials with Stackel po- 
tentials are required (see for instance de Bmyne et al 12000b to 
obtain more reliable estimates of the vertical tilt angle. 
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